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As part of a big national scientific network ’Pathways to Habitable Worlds’ the formation of 
planets and the delivery of water onto these planets is a key question because water is essential for 
the development of life on a planet. In the first part of the paper we summarize the state of the art 
of planet formation - which is still under debate in the astronomical community - before we show 
our results on the formation of planets. The outcome of our numerical simulations depends a lot on 
the choice of the initial distribution of planetesimals and planetary embryos after the disappearance 
of gas in the protoplanetary disk. In addition we take into account that some of these planetesimals 
of sizes in the order of the mass of the Moon already contained water; the quantity depends on the 
distance to the Sun - too close and the bodies are dry, but starting from a distance of about 2 AU they 
can contain substantial amounts of water. Our assumption is that the gas giants and the terrestrial 
planets are already formed when we check the collisions of the small bodies containing water (in 
the order of a few percent) with the terrestrial planets. We thus are able to give an estimate of the 
respective contribution to the actual water content (of some Earth-oceans) in the mantle, in the 
crust and on the surface of Earth. In the second part we discuss how the formation of larger bodies 
after a collision may happen in detail, because the outcome depends on different parameters like 
the collision velocity, the angle of the impact and the materials involved. The detailed results were 
accomplished with the aid of SPH (Smoothed Particle Hydrodynamics) simulations, a sophisticated 
numerical method to simulate these events. We briefly describe this method and show different 
scenarios with respect to the formed bodies, possible fragmentation and the water content before 
and after the collision. 

In an appendix we discuss the detection methods for extrasolar planets around other stars for which 
the current number of known ones is already close to 2000. 


STATE OF THE ART 

The formation of terrestrial planets is an outstanding 
problem in astronomy which has - up to now - not been 
solved satisfactory. Especially difficult is to form a low- 
mass Mars at about 1.5 AU and a very dense Mercury 
inside the orbit of Venus. All of the theories in this re¬ 
spect assume that the gas giants in the outer Solar Sys¬ 
tem were formed first from a gas dominated protoplan¬ 
etary disk in relatively short time [1] and that the inner 
planets were formed later. Globally one can say that 
the outcome of every simulation is highly dependent on 
the initial conditions for the distribution of mass in the 
early Solar System. The Mars problem has led to differ¬ 
ent approaches concerning the scenarios for the growth of 
planetesimals to planetary embryos to protoplanets and 
finally to planets. Because of the small mass - 1/10 of 
Earth’s mass for Mars (and Mercury) - they are regarded 
as not fully developed planets, but rather as protoplanets 
which did not accumulate mass in later stages of the early 
history of the Solar System. While for Earth (and also for 
Venus) a formation time of about 100 million years seems 


appropriate, Mars should have formed within 10 million 
years only. One recent theory is the so-called GrandTack 
model [2]: 

A disk (up to 3 AU) of relatively dry planetesi¬ 
mals in low-eccentricity orbits resides inside an al¬ 
ready fully formed Jupiter at about 3.5 AU and Sat¬ 
urn at about 5 AU still accumulating material. The 
inward migration of the gas giants caused by the still 
present gas drag truncates the original disk to an 
outer edge of 1 AU. When Jupiter was at 1.5 AU 
and Saturn (now fully formed) was caught in a 2:3 
mean motion resonance at about 2 AU the migration 
was reversed. Without a deeper explanation in the 
original paper it was stated that ’ ... according to our 
results Jupiter tacked at 1.5 AU and then migrated 
outward owing to the presence of Saturn’ [5]. During 
this stage of migration outward, outer wet material 
was scattered to the region inside Jupiter’s orbit, so 
that finally the distribution of the asteroid main belt, 
composed of a dry inner region (S-type asteroids), and 
a wet outer region (C-type asteroids), as we see it to¬ 
day, was reached, with the gas giants more or less in 
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their present positions. Now the formation of the ter¬ 
restrial planets to their final masses and compositions 
took place in this inner region, which already started 
during the outward migration of the gas giants. 

Prior to that interesting but very special scenario dis¬ 
cussed before, where many assumptions needed to be 
made, a simple simulation [3], starting with 400 equally 
massive bodies (0.005 Earth-masses each), randomly dis¬ 
tributed in an annulus between 0.7 and 1 AU with low 
eccentricities and inclinations, has shown a surprising re¬ 
sult. The distribution of the terrestrial planets after long 
term integrations turned out to be the one which we know 
from our present-day Solar System: a small Mercury in¬ 
side, and a small Mars outside the orbits of the ’twin 
planets’ Earth and Venus close to their present positions. 

A recently published theory [3] uses these results and 
explains the existence of a low-mass Mars by straight 
forward integrations with a special assumption for the 
distribution of about 150 planetary embryos and about 
1000 planetesimals inside a disk extending from 0.5 AU 
to 4 AU in the presence of Jupiter and Saturn in their 
current positions. Their major point is that they assume 
a distribution of planetesimals and embryos being much 
lower around 1.5 AU, where Mars can be found today. 
According to an article by Jin et al. [5] this local min¬ 
imum can be explained by the evolution of the surface 
density distribution of the solar nebula during the first 
million years. But as pointed out by Raymond et al. 
[H] and Raymond et al. [7] this dip is too narrow to cut 
off Mars’ accretion. Given the broadness of the lower- 
density gap used by Izidoro et al. [4] the results shown 
there should be taken with caution. 

It should be emphasized that for the last 20 years the 
interest in studies about the formation of planets and 
especially terrestrial planets with respect to their habit¬ 
ability is growing quasi exponentially and the literature 
of only the last years is huge. There are some summa¬ 
rizing review articles on this topic which more or less 
discuss the latest progress, but all of these papers are 
very much directed towards the authors’ own research 
and leave many questions out. Just a few of them should 
be mentioned: Agnor and Asphaug [8], Raymond et al. 
[6] , Marcus et al. [9] , Agnor and Lin [10] , Morbidelli et al. 
[TT] , Raymond et al. |7| , Lineweaver and Chopra [T^l and 
Giidel et al. [Bj. 


OUR METHODS 

One main constraint in all the cited simulations is the 
number of integrated planetesimals and planetary em¬ 
bryos which was always limited according to the integra¬ 
tion times on single CPUs. Although some interesting 
codes are available, where up to about 1000 bodies can 
be integrated in the Newtonian framework (e.g. Cham- 



FIG. 1. Two colliding bodies at the moment of the first en¬ 
counter: the two important parameters are the impact veloc¬ 
ity and the impact angle a besides their radii Rp and Rp. 

bers M), these simulations are not precise in a very im¬ 
portant fact: the collisions of two bodies which merge are 
not precisely computed, which causes a major drawback. 
In this connection we were and are still working into dif¬ 
ferent directions: first of all we avoid simplifications for 
collisions and the merging of bodies insofar as that we 
compute the close encounters up to the moment when 
the two bodies involved are actually touching; secondly 
we extended the programs to be able to integrate several 
thousands of interacting bodies by writing codes running 
on graphic cards. In the following we will concentrate on 
the first point. 


The Encounters and the Collisions 

As mentioned one crucial step in the process of forming 
planets are the collisions of the early Solar System bodies 
and their possible merging. Ghecking such an encounter 
in detail one can see that as well as the impact velocity, 
the impact angle is a different one when we assume a 
collision happening already at a larger distance than the 
real physical one as depicted in Fig.l and Eig.2. 

The increasing velocity visible in Figj^ (green line) 
leads to a sharp cusp at the closest distance (red line). In 
fact the collision between two extended bodies (e.g. two 
Geres-sized bodies) and the following merging happens 
well before the moment of entering the innermost small 
circle (green, ro,i for two bodies of the mass of I/IO of 
Geres) in Figj^and depends on the diameters of the two 
colliding bodies. The second orbit in this figure seems 
to be tangential to the blue circle representing the col¬ 
lision distance (rgy, depicted in Figfor two arbitrary 
objects). 
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3 examples of flybys (collisions) 



FIG. 2. Planetocentric encounters with three different veloc¬ 
ities and impact angles when entering a distance of 0.001 AU 
(red circle): a high velocity encounter (magenta), a moder¬ 
ate velocity encounter (blue) and a relatively far away fly by 
(yellow) are depicted (for more details see the text). 


It has to be mentioned that different velocities and col¬ 
lision angles lead to different outcomes, e.g. disruption, 
partial accretion and smooth merging, depicted in Fig|^ 
Seven different masses were chosen to study the colli¬ 
sion velocities, ranging from 1/10 of Ceres to 10 times 
the mass of the Moon (for explanation see table 1); note 
that very high velocities only occur in the 1/10 Ceres sce¬ 
nario. The lines in the illustration refer to the boundaries 
of different collision outcomes as given in the paper by 
Leinhardt and Stewart m- erosion of the target above 
the dot-dashed curve, below it partial accretion to the 
left and hit-and-run scenarios to the right of the thick 
vertical line. Finally a thin region of perfect merging can 
be found just above the lowermost dashed line, while be¬ 
low it (where v/vesc < 1) no collision is possible at all. 
A more qualitative illustration of this behaviour can be 
found in Fig|^ 

It is evident from Figthat the assumption of a sim¬ 
ple inelastic merging for all collisions is not true and an 
oversimplification. Except for some notable exceptions 
(e.g. Genda et al. US]) nobody has yet accounted for 
this rather diverse range of collision outcomes. We will 
concentrate on this issue in the subsequent sections when 
we discuss results of SPH simulations of collisions. Com¬ 
monly the scenario of complete merging is used as an 
approach to the problem; in the next section we show 
results of our own computations where we made these 
same ad hoc assumptions. 
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FIG. 3. Planetocentric distances (red) and collision velocities 
(green) plotted versus time. A high velocity encounter (upper 
graph) and a low velocity encounter (lower graph) are shown 
(for more details see text). 
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TABLE I. N-body simulation scenarios - m is the planetes- 
imal mass, i?imp (in 10® m) denotes the mutual distance of 
the planetesimals’ barycenters upon impact for rniii the Hill 
radius at lAU and zero eccentricity, and Uesc the target’s 
surface escape velocity (see text); All scenarios include 750 
planetesimals and are integrated for 1 Myr. 

RESULTS OF NEW COMPUTATIONS 

As already mentioned most current scenarios of planet 
formation are lacking two major points: the ’correct’ ini¬ 
tial distribution of the planetesimals/planetary embryos 
and realistic collision outcomes. Whereas the second 
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Impact angle [deg] 
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FIG. 4. Collision velocities (normalized by the surface escape velocities) versus impact angle for the collision of differently sized 
bodies from mceres/10 to 10 rriMoon (see detailed explanation in table 1 and in the text). 


point will be taken into account in future models only, 
the first point was considered by using the outcome of the 
Grand Tack model: we distributed 100 dry planetesimals 
and planetary embryos randomly (with small inclinations 
and eccentricities) in the range of 0.5 AU < a < 2.2 AU. 
Jupiter and Saturn were in their current positions, while 
Uranus and Neptune were not present in the simulated 
debris disk. Whenever two small bodies suffer from a col¬ 
lision the newly formed body replaces them and a new 
smaller body (with low water content) is added from the 
outer ’wet’ belt of planetesimals. This - strictly speak¬ 
ing - somehow artificial procedure was adopted due to 
numerical issues and aims to simulate the water trans¬ 
port to the inner dry region. The reason is to keep the 
number of bodies limited to 100, which is in our numer¬ 
ical approach, the Lie-series technique HI], a number of 
bodies that allows to stay in a reasonable amount of com¬ 
puting time. Without going into the details of the nu¬ 
merical integration procedure it should be noted that the 
automatic step-size control makes the computation more 
accurate compared to other methods [H] (especially con¬ 
cerning close encounters). 

In the respective Figj^we show the summarized results 
of 16 independent simulations with different initial condi¬ 
tions with Jupiter as sole perturbing planet of the cloud 
of planetesimals (embryos). The small red points mark 
the initial conditions for the 100 primary bodies, the dots 
represent the binned outcome after 1.5 Myrs. Two dif¬ 
ferent interpolations show the peak at a=l AU which is 
the position of Earth; the body with the largest mass has 
approximately half the mass of present-day Earth. 

Another run was performed with only six independent 
simulations, shown in Figj^ but with Jupiter and Sat¬ 
urn as perturbers of the belt of bodies inside Jupiter. 
Here one can observe two accumulation maxima, repre¬ 
senting protoplanets, due to repeated merging processes: 


one at a=l AU and another one with a=1.5 AU. The 
second peak with a considerable lower mass (about 1/10 
of Earth’s mass) is exactly at the location of Mars. 

One representative example out of the 22 different sim¬ 
ulations is depicted in Fig[^ After 1.5 million years a 
planet of ca. 1/2 of Earth’s mass is formed in a dis¬ 
tance of almost 1 AU (planet A). This planet has a 
low eccentricity and a very small inclination. Four more 
smaller planets (embryos with a mass around 10 times 
the Moon’s mass) are formed: 

1. Planet B: at 0.59 AU, a Mercury like object with a 
relatively large eccentricity but a small inclination 

2. Planet C: at 0.98 AU, an object very close to planet 
A with almost the same eccentricity but a slightly 
larger inclination 

3. Planet D: at 1.5 AU another embryo with eccen¬ 
tricity and inclination comparable to planet C in 
the distance of Mars in our Solar System 

4. Planet E: at 1.7 AU, an object with an eccentricity 
comparable to D but with a large inclination. 

The planets A and C were separately integrated for 
another 10 Myrs and it turned out that after 6 Myrs 
they merge to a bigger body at approximately 1 AU. This 
shows that to form Earth-sized planets the time scale is 
at least in the order of tenth of Myrs. 

The Watertransport 

As mentioned before the inner part of the primordial 
disk of planetesimals is dry so that we need to find an 
explanation for the significant amount of water on Earth 
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FIG. 5. Synopsis of simulations after 1.5 Myrs only with Jupiter for 16 different initial conditions distributed according to the 
outcome of the GT scenario (red points). The plot shows the merging planetesimals’ semimajor-axes versus the mass in solar 
masses. A peak at 1 AU, corresponding to the position of Earth, is clearly visible. 



FIG. 6. Synopsis of simulations after 1.5 Myrs with Jupiter and Saturn for 6 different initial conditions distributed according 
to the outcome of the GT scenario. The plot shows the merging planetesimals’ semimajor-axes versus the mass in solar masses. 
Two peaks are visible, one at 1 AU and another one at about 1.5 AU. These two peaks correspond approximately to the 
positions of Earth and Mars. The magenta line marks a mass of 10 times the Moon’s. 
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FIG. 7. Results of planet formation after 1.5 Myrs: semimajor 
axis versus eccentricity (upper picture) and versus inclination 
(lower picture) respectively. The size of the points is propor¬ 
tional to their mass; the largest body formed at 1 AU and has 
about 30 times the mass of the Moon. The color code indi¬ 
cates the mass of water on the formed planets in Earth-oceans 
(for detailed explanation see text). 



FIG. 8. Image of two colliding bodies of the size of Geres 
computed with our SPH-code. The target has no water at all, 
the projectile has a layer of water on the surface (30 percent 
of its total mass). 


inside the snow line. This border is located at approx¬ 
imately 2 AU and it is by definition the minimum dis¬ 
tance to the Sun for which water can be present in the 
liquid state and also as ice. The collision behavior (see 
Figj^for an illustrative example) of small bodies coming 
from outside the snow line regarding their water con¬ 
tent is crucial to answer this question. They have to be 
modeled with specially designed programs like our SPH 
(Smoothed Particle Hydrodynamics) code, which will be 
discussed in detail in the next chapter. In these simula¬ 
tions we fully take into account the complicated process 
of (giant) collisions. For our N-body integrations on the 
other hand we just assume perfect inelastic merging of 
two colliding bodies and determine the respective water 
loss at this stage. The outcome of one such scenario as an 
example is shown in Fig[^ We can see that planet A with 
a mass of ~ ^/‘^MEarth has a water content of at least 
two oceans, which is in accordance with the estimated 
4-8 oceans (surface and mantle combined) of present- 
day Earth. Rather surprising is the water content of 
the innermost planet B, located at Mercury’s position 
but, we did not take into account the close proximity 
to the central star, which may have made Mercury dry. 
Many questions are still open and we are still far away 
from understanding the full complexity of planet forma¬ 
tion, especially in connection to planet Earth, the ’BLUE 
DOT’ in our Solar System. 

SIMULATING COLLISIONS WITH SPH 

It is well established that planet formation results from 
a sequence of collision events between protoplanetary 
bodies of different sizes over an extended period of time. 
This collisional growth of planetesimals is investigated 
by means of dynamical evolution studies that are based 
on N-body simulations. The underlying collision model 
that leads to growth is typically simplified and based on 
conserving linear momentum in perfect inelastic merging 
[e.g., UMI] or deploying simple fragmentation models 
| 22 ) . As this approach neglects potential fragmentation 
in partial accretion and hit-and-run encounters, it misses 
the actual physics of collisions [221 HI]- Rather, the re¬ 
sulting configuration after a collision depends on the in¬ 
volved bodies’ mass ratios, encounter velocities, and im¬ 
pact angles. It can be categorized in one of the four 
regimes efficient accretion/perfect merging, partial accre¬ 
tion, hit-and-run, and erosion and disruption |24j . Ag- 
nor and Asphaug for example, study giant collisions’ 
accretion efficiency by simulating collisions using SPH 
(projectile and target masses of 0.1 Mq); Marcus et al. 
[26] extend disruption criteria for giant impacts up to 
masses of 10 M 0 . As a consequence, one might expect 
consistently overestimated masses of the surviving bod¬ 
ies as supported e.g., by Alexander and Agnor [15] who 
investigated different simplified merging models and find 
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the same number of big bodies on similar orbits, but de¬ 
pending on the merging assumption the bodies’ masses 
increase towards higher merging rates. According to the 
computing power available at that time, Alexander and 
Agnor simulate in 2D for up to 10'^ yr and deploy a 
very rudimentary fragmentation model which, however, 
predicts up to four surviving fragments and a remaining 
core. Kokubo and Genda [27] aim for a merging criterion 
and derive an empirical formula for the critical impact 
velocity that discriminates merging and hit-and-run col¬ 
lisions. They base their collision data on SPH collision 
simulations (resolution 20k-60k SPH particles). In sim¬ 
ulations of protoplanetary systems for 3 x 10® yr they 
find that only about half of the collisions lead to accre¬ 
tion. Nevertheless, no significant influence of the merging 
model on planetary growth timescales and the masses of 
the largest planets is found. 

Agnor et al. |23j analyzed the bodies’ angular momenta 
and show that assuming perfectly inelastic merging can¬ 
not be sustained because it would lead to rotationally 
unstable bodies. Spin rates resulting from the criterion 
by Kokubo and Genda m are about 30 % lower than in 
the perfect merging case. 

All in all, 40-50% of giant collisions are found to result 
in a non-merged configuration [251 EH EH]. 

As we are interested not only in the formation of plan¬ 
ets, but also in water transport we decided to perform 
collision simulations including water content on the col¬ 
liding bodies and consider effects of different material 
strengths of rocky and icy components. 


Physical model 


Most giant impact studies [e.g., ES] use strengthless 
material based on self-gravity dominating the material’s 
tensile strength beyond a certain object size [400 m in ra¬ 
dius, EOj ■ lu Maindl et al. EH we investigate collisions of 
objects close to this size limit and compare strengthless 
material simulation (“hydro model”) to modeling ma¬ 
terial strength with associated elasto-plastic effects and 
damage/brittle failure (“solid model”) which we will out¬ 
line in this chapter. 

We use the Tillotson |HH| equation of state (EOS) as 
formulated in Melosh |33| for modeling the bodies’ ma¬ 
terials. It depends on 10 material constants po, A, B, a, 
6, a, P, Eq, Eiy, and E^v and distinguishes two domains 
depending upon the material specific energy E. In com¬ 
pressed regions with density p > po) and E lower than 
the energy of incipient vaporization E{y the EOS reads 


P = 


1 + E/(Eo772)J 


pE -\- A/ix ~\~ Bp.^ (1) 


the EOS is given by 


P = apE+ 


hpE 

1 + E/iEop^) 


Apr 

e/3(po/p-i) 


g-a(po/p-l)^ 


( 2 ) 


Linear interpolation between the pressures obtained 
via Q and ([^ gives the P in the partial vaporization 
regime Ejv < E < Ecv (refer to Melosh (33] for a more 
detailed description). 

Gonservation of mass, momentum and energy describe 
the dynamics of a solid body according to the theory of 
continuum mechanics, cf. Schafer et al. |34j . The conti¬ 
nuity equation provides conservation of mass and reads 
in Lagrangian form (Einstein notation) 


dt 


dv° 

dx'^ 


= 0 . 


Momentum is conserved by 


du“ 

dt 


1 

p dx^ 


= -P6°‘^ + S°‘^ 


with the stress tensor (t“^ that can be split between 
the pressure P along its diagonal and and the deviatoric 
stress tensor S°‘^ (5“^ denotes the Kronecker delta). En¬ 
ergy conservation reads 


^ ^ _P^ ^ 

dt p dx°‘ p 

where denotes the strain rate tensor as given in 
While the EOS provides an analytic function for the 
pressure, the time evolution of the deviatoric stress tensor 
S°‘^ is needed to completely describe the dynamics of a 
solid body. We use Hooke’s law to define the constitutive 
equation as 


dS°‘^ 

dt 



1 

3 


gapp-y 




where p represents the shear modulus and the strain 
rate tensor. The last two terms depend on the rota¬ 
tion rate tensor and guarantee that the constitutive 
equations are independent from the material frame of ref¬ 
erence. and are given by 


pap 


2\dxP dx°‘J ' 


1 /'dv°‘ 

2 



■ (3) 


Equations 0 to ^ describe the dynamics of an elastic 
solid body. We follow the approach by Benz and Asphaug 
[35] to model plastic behavior and use the von Mises yield 
criterion where the material yield stress Y limits the de¬ 
viatoric stress. Hence, a transformed deviatoric stress 
introduced: 


qaP 

‘^vM 


= min 


2K" ■ 

3 5'a/35'a/3’ 




where ry = p/po and pT = rj — 1. For expanded states 
{E greater than the energy of complete vaporization Ecv) 















We model tensile failure (fracture) by implementing 
a damage model based on the fragmentation model of 
Grady and Kipp |36j and the SPH implementation by 
Benz and Asphaug [35]. The model is based on assigning 
flaws to SPH particles that are activated once the strain 
exceeds a certain level e. Activated flaws convert into 
cracks. The ratio of cracks to flaws defines the damage 
value d £ [0,1] assigned to each SPH particle. The devi- 
atoric stress and tension decrease proportionally to 1 — d 
and vanish for completely damaged material (d = 1). Ini¬ 
tially, the number n of flaws that can be activated by a 
strain level of e are distributed according to a Weibull dis¬ 
tribution with material parameters k and m: n(e) = fc e"* 

m- 

Table [H] gives the Tillotson EOS parameters, the shear 
modulus /X and yield stress Y which we use in our simu¬ 
lations. The values are adopted from Benz and Asphaug 
[88j following their reasoning on approximating po: set¬ 
ting A equal to the bulk modulus, and B — A. 

We use Weibull flaw distribution parameters of basalt 
based on directly measured data by Nakamura et al. |89] 
and for ice we adopt the values mentioned in Lange et al. 
[dOj. rUbasalt ~ 16, A^basalt ~ 16 7 ^ice ~ 9.1, kicQ = 

10^6 m-3. 


Fragmentation and water loss 

The presented study was performed with our own 
smoothed particle hydrodynamics (SPH) code which is 
based on the physics described in the previous section 
with the addition of self-gravity and a tensorial correction 
ensuring first-order consistency as described in Schafer 
et al. [33] . The SPH code has been introduced in Maindl 
et al. [31], Schafer mi, Maindl et al. [33] and used for 
simulating high-resolution surface impacts [331 El] and 
planetesimal-to-protoplanet collisions [SiiSiEg. 

Our scenarios presented here represent two colliding 
objects that have a mass of Mceres each. The target is 
assumed to consist of a rocky basalt core (70mass-%) 
and a mantle/shell of water ice (30mass-%) whereas the 
projectile is purely rocky (basalt). In total we simulated 
42 scenarios parameterized by collision velocities vq be¬ 
tween 0.95 and 5.88 two-body escape velocities Uesc and 
impact angles a between 0° (head-on) and a flyby. As 
we are interested in an overall fragmentation and water 
transport pattern we use a relatively low resolution of 
20k SPH particles. 


Collision outcome 

Due to the typically high degree of fragmentation after 
impacts in the solid model and in order to get statisti¬ 
cally relevant physical fragment properties we limit our 
study to fragments with masses Wfrag corresponding to at 


least 20 basalt SPH particles. Given the scenarios’ single 
basalt SPH particle mass of 1.2 x 10^^ kg this results in 
the condition rufrag > 2.4 x 10^® kg while the remaining 
mass will be considered as lost debris. This limits the res¬ 
olution of our results. Table [III] states the parameters of 
the 42 different collision configurations (impact velocity 
and angle) simulated in Maindl et al. [3T] along with the 
resulting numbers Wrag of such “significant fragments” 
after the collision. Each a-Wrag column pair corresponds 
to one initial impact parameter. The bodies are started 
five mean diameters apart to allow the particle distribu¬ 
tion to settle due to mutual tidal forces. Because of mu¬ 
tual gravitational interaction different encounter speeds 
result in different orbital deflection and hence different 
actual impact angles a. In case of the largest initial im¬ 
pact parameter (rightmost column in Tab. HI) only the 
slowest encounter velocity leads to an impact; otherwise 
the scenario results in a flyby and hence cannot be as¬ 
signed an impact angle. The latter show two surviving 
fragments which are exactly the original bodies. 

Fig. [| relates our fragmentation results to an ana¬ 
lytic collision outcome model for strengthless planets pre¬ 
sented by Leinhardt and Stewart [T3]. Overall, the fV^ag 
numbers correspond to the model, with some exceptions 
at the boundaries of the collision outcome regimes. They 
suggest that the onset of the hit-and-run regime actu¬ 
ally happens at a higher impact angle than predicted by 
the strengthless model. At higher angles the fVfrag = 6 
bubble actually represents a hit-and-run with four very 
small fragments which is in accordance to the analyti¬ 
cal model, whereas the fVfrag = 7 outcome is actually a 
merging event with 6 very small fragments and suggests 
the merging limit to be slightly too low in the analytic 
model. Please refer to Maindl et al. [311 Sg for a more 
thorough discussion of the detailed outcomes. 


Water content 

Tracking overall water content retained in significant 
fragments after the collision and comparing to the water 
content on the original bodies gives an estimate on the 
loss of volatiles such as water ice. Figure [^ shows how 
much of the water is still on significant fragments after 
the collision. The remaining water is on the debris and 
lost in our model. For a < 20° and uq ^ l-Suesc most of 
the water ice stays on the survivor, in strongly inclined 
hit-and-run scenarios more than 80 % of the water stays 
as well. In general more water is lost for “more violent” 
impacts—an increasing amount of water ice is lost in de¬ 
bris for small collision angles and high velocities with a 
more significant dependency on velocity than on the an¬ 
gle. 

Notable features in Fig. [^include downward spikes oc¬ 
curring in the low velocity scenarios (vq < 1.36Uesc)- A 
more thorough investigation reveals that these are caused 
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TABLE II. Tillotson EOS parameters, vaporization energy levels, shear modnlns fi, and yield stress Y, cf. Benz and Asphang 
[38]. A — B is set equal to the bulk modulus. 



Po 

(kg/m^) 

A 

(GPa) 

B 

(GPa) 

Eo 

(MJ/kg) 

Eiv 

(MJ/kg) 

Ecv 

(MJ/kg) 

a 

b a 

P 

P 

(GPa) 

Y 

(GPa) 

Basalt 

2700 

26.7 

26.7 

487 

4.72 

18.2 

0.5 

1.50 

5.0 

5.0 

22.7 

3.5 

Ice 

917 

9.47 

9.47 

10 

0.773 

3.04 

0.3 

0.1 

10.0 

5.0 

2.8 

1 


TABLE III. Number of significant surviving fragments A^ag after 2000 min simulation time. The scenario parameters 
collision velocity vo given in units of Vssc (two-body escape velocity) and the impact angle a. The wo (averaged) and a 
are taken from Maindl et al. [STji 8- missing impact angle indicates a flyby scenario (see text). 

are the 
values 

^0 

[t^esc] 

a 

n 

-^frag 

a 

ri 

A^frag 

a 

n 

-^frag 

a 

n 

A^frag 

a 

n 

-^frag 

a 

n 

A^frag 

0.95 

0 

1 

14 

1 

23 

1 

25 

1 

31 

1 

62 

7 

1.32 

0 

1 

11 

1 

21 

1 

40 

1 

48 

2 

- 

2 

1.36 

0 

1 

11 

1 

20 

1 

40 

2 

48 

2 

- 

2 

2.12 

0 

1 

12 

1 

25 

5 

50 

2 

62 

6 

- 

2 

3.04 

0 

39 

12 

52 

25 

22 

53 

2 

67 

2 

- 

2 

3.97 

0 

41 

13 

67 

27 

35 

55 

2 

72 

2 

- 

2 

5.88 

0 

45 

13 

60 

28 

61 

58 

2 

72 

2 

- 

2 


by a single surviving body spinning rapidly losing debris 
including large portions of its surface water ice. A more 
detailed discussion of fragment counts and individual sce¬ 
nario descriptions is available in Maindl et al. ED- 


Future work 

Based on earlier results of simulating small to mid 
scale collisions of planetesimals at moderate energies 
EH He] we presented first conclusions on fragmentation 
and volatile loss. Rather than analyzing total fragment 
counts we identified surviving “significant fragments” 
and their volatile content. While the collision outcomes 
qualitatively agree with analytic giant collision outcome 
models of strengthless planets [as given e.g., in US], 
we observe some minor shifts of boundaries between the 
merging, erosion, and hit-and-run regimes. We further 
observe significant water loss—up to 60%—for faster 
and/or less inclined impacts. This might have a major 
impact on future planet formation simulations as most 
of them currently rely on a 100 % water transfer when¬ 
ever a collision occurs. Incorporating a realistic model 
of volatile transport and loss may result in reduced wa¬ 
ter abundances of formed terrestrial planets by a factor 
of 5-10 H?]- Envisioned future studies will include frag¬ 
ment dynamics determining their ability to escape the 
system’s Hill sphere connected with working towards a 
comprehensive model for the fate of volatiles such as wa¬ 
ter ice and including it in N-body simulations for planet 
formation. 


Appendix: The Discovery of Exoplanets 

Up to now we have knowledge of 1164 planetary 
systems, where 473 of them host multiple compan¬ 
ions, giving a total of 1855 known extrasolar planets 
{www.exoplanet.eu, last updated 2014). The believe in 
other worlds similar to our own dates back to the ancient 
Greeks like Democritus (~ 460-370 B.C.) and Epicurus 
(341-270 B.C.). 2000 years later, in the sixteenth cen¬ 
tury, the Italian philosopher Giordano Bruno, who was a 
supporter of the Gopernican theory, put forward the idea 
that the stars are similar to the Sun and could also be 
accompanied by planets. 

It is interesting to note that the first confirmed plan¬ 
ets, orbiting a pulsar |48j . were rather small compared to 
most later discoveries. Only 5 years later the first planet 
around a main sequence star, a hot Jupiter with a 4.2-day 
orbit (51 Pegasi), was discovered by Mayor and Queloz 

En¬ 

in the following we summarize the different techniques 

for discovering extrasolar planets: 

• radial velocity technique (RV) 

• transit photometry (TP) 

• microlensing (ML) 

• direct imaging (DI) 

The first three are indirect methods, making it impos¬ 
sible to directly separate the orbiting planet from the 
host star. Eurthermore one has to distinguish ground 
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Impact angle a [deg] 

0 15 30 45 60 90 



0 -'-^^^i^^^i^^-'-'-'-'-^^^^- lJ 

0 0.2 0.4 0.6 0.8 1 

Impact parameter b = .sin(a) 

FIG. 9. Number of significant surviving fragments for the scenarios’ collision velocities Vo given in units of Ueac (two-body 
escape velocity) and impact angles a. Different areas in the graph correspond to erosion (red, shaded), partial accretion (white, 
a < 30°), hit-and-run (white, a > 30°), and perfect merging (blue, hashed) as given by the analytic model for strengthless 
planet collisions in Leinhardt and Stewart Fig. IIA. See text for discussion. 



FIG. 10. Gumulated water retained in significant fragments after the collision. A value of one means all water in the system 
is on the surviving significant fragments. The different curves correspond to different collision velocities vq in units of the 
two-body escape velocity Vbsc- See text for discussion. 


based detections and observations via space telescopes, a 
method up to now only applied to the TP technique. 

1. RV This method, also called Doppler technique, 
uses the Doppler effect for analyzing the motion 
and properties of the star and planet. Because 
both, the planet and the star, are orbiting their 
common center of mass, spectral lines are shifted 
versus shorter wavelengths (blue) when the host 
star moves towards the observer. When the planet 
moves towards the observer on the contrary, the 
star is moving away, and consequently the spec¬ 
tral lines of the star are shifted versus longer wave¬ 
lengths (red). This change happens periodically 
since the planet orbits the star with a certain pe¬ 
riod. This allows us to obtain a so-called RV curve, 
where one can very well determine the period and 
even the orbital eccentricity from. 

2. TP Like Sun, Moon and Earth are aligned during 


a solar eclipse, observer, planet and its host star 
are in one straight line during a TP event. The ex¬ 
trasolar planet is moving in front of the stellar disk 
just like we can observe it on Earth for Mercury 
and Venus every now and then. Knowing the real 
diameter of the star via spectroscopic observations 
and also in addition its RV curve (typically from 
follow-up observations) one can determine the pe¬ 
riod of the orbit, the planet’s diameter, the orbital 
eccentricity and in some instances also the inclina¬ 
tion. Having obtained these values it is straight 
forward to estimate the exoplanet’s mass and bulk 
density with the aid of the 3'"'^ law of Kepler. These 
parameters are important because at least for most 
examples they tell us whether it is rather a gaseous 
planet or a rocky one. 

3. ML When a distant star is almost aligned with a 
massive compact foreground star hosting a planet, 
the distortion of light due to its gravitational field 
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may lead to two unresolved images and thus to an 
observable magnification. This transient brighten¬ 
ing depends on the mass of the object which acts 
as a lens. It is important to mention that this is al¬ 
ways a unique event because the almost alignments 
of the observer on Earth and the two other objects 
are quite rare events. In addition the time of pos¬ 
sible observations depends on the relative proper 
motion between the background ’source’ and the 
foreground ’lens’ object. 

4. DI Because of the big difference in luminosity be¬ 
tween the star and a planet this technique has so 
far been successful only for big planets relatively 
far away from their host star. 

For the last ten years photometric observations from 
space are carried out as well: the first larger instrument 
was the CoRoT (COnvection, Rotation and planetary 
Transits) satellite, started in 2006 by ESA and was work¬ 
ing well for about 8 years. With a 27cm telescope on 
board thousands of stars were investigated, leading to 
hundreds of planetary candidates and about 40 confirmed 
exoplanets. NASA’s KEPLER satellite observes with a 
95cm telescope and was launched in 2009. It is still work¬ 
ing and the bigger diameter of the telescope allowed for 
observations of hundreds of thousands of stars, resulting 
in more than 1000 planetary candidates and several hun¬ 
dreds of confirmed planets by later measurements from 
the ground. When speaking about candidates and con¬ 
firmed planets it should be mentioned that observations 
from space alone are not enough, ground based observa¬ 
tions, mainly spectroscopic ones for obtaining additional 
RV curves, need to be done. With the TP method one 
may be observing the dimming of light caused by star 
spots, an apparently close by but indeed far away back¬ 
ground star or even a close binary, leading to a false¬ 
positive planet detection. 

To summarize one can say that the last twenty years of 
research on extrasolar planets opened a new and exciting 
chapter in the field of astronomy, and also for the human 
civilization as a whole, trying to answer the big question: 
ARE WE ALONE? 
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